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Page 22, mid-page, correct the following code: 
Old code: A = 6.4837D0 
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SPECIAL NOTES 


1. API PUBLICATIONS NECESSARILY ADDRESS PROBLEMS OF A GENERAL 
NATURE. WITH RESPECT TO PARTICULAR CIRCUMSTANCES, LOCAL, STATE, 
AND FEDERAL LAWS AND REGULATIONS SHOULD BE REVIEWED. 


2. APL IS NOT UNDERTAKING TO MEET THE DUTIES OF EMPLOYERS, MANU- 
FACTURERS, OR SUPPLIERS TO WARN AND PROPERLY TRAIN AND EQUIP 
THEIR EMPLOYEES, AND OTHERS EXPOSED, CONCERNING HEALTH AND 
SAFETY RISKS AND PRECAUTIONS, NOR UNDERTAKING THEIR OBLIGATIONS 
UNDER LOCAL, STATE OR FEDERAL LAWS. 


3. INFORMATION CONCERNING SAFETY AND HEALTH RISKS AND PROPER 
PRECAUTIONS WITH RESPECT TO PARTICULAR MATERIALS AND CONDI- 
TIONS SHOULD BE OBTAINED FROM THE EMPLOYER, THE MANUFACTURER 
OR SUPPLIER OF THAT MATERIAL, OR THE MATERIAL SAFETY DATA SHEET. 


4. NOTHING CONTAINED IN ANY API PUBLICATION IS TO BE CONSTRUED AS 
GRANTING ANY RIGHT, BY IMPLICATION OR OTHERWISE, FOR THE MANU- 
FACTURE, SALE OR USE OF ANY METHOD, APPARATUS, OR PRODUCT 
COVERED BY LETTERS PATENT. NEITHER SHOULD ANYTHING CONTAINED 
IN THE PUBLICATION BE CONSTRUED AS INSURING ANYONE AGAINST 
LIABILITY FOR INFRINGEMENT OF LETTERS PATENT. 


5. GENERALLY, API STANDARDS ARE REVIEWED AND REVISED, REAF- 
FIRMED OR WITHDRAWN AT LEAST EVERY FIVE YEARS. SOMETIMES A ONE 
TIME EXTENSION OF UP TO TWO YEARS WILL BE ADDED TO THIS REVIEW 
CYCLE. THIS PUBLICATION WILL NO LONGER BE IN EFFECT FIVE YEARS 
AFTER ITS PUBLICATION DATE AS AN OPERATIVE API STANDARD OR, WHERE 
AN EXTENSION HAS BEEN GRANTED, UPON REPUBLICATION. STATUS OF THE 
PUBLICATION CAN BE ASCERTAINED FROM THE API PUBLICATIONS DEPART- 
MENT [TELEPHONE (202) 682-8000]. A CATALOG OF API PUBLICATIONS AND 
MATERIALS IS PUBLISHED ANNUALLY AND UPDATED QUARTERLY BY API, 
1220 L STREET, N.W., WASHINGTON, D.C. 20005. 
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FOREWORD 


A simplified correlation for the equilibrium bubble point pressure, vapor pressure, of 
commercial natural gas liquids, NGLs, has been developed. This correlation was developed 
in response to the industry’s need for a simplified correlation to be used when adjusting 
densities from flowing conditions to a set of standard conditions. Methods for improving 
the accuracy and extending the correlation are also discussed. 

The purpose of this procedure is to provide a simplified means of estimating flowing 
equilibrium vapor pressures of various Natural Gas Liquids, using relative density (60°F/ 
60°F) and flowing temperature inputs. This procedure is proposed for use on GPA Standard 
Specification 2140, commercial HD-5 Propane, Iso-Butane, Normal Butane, and 12-14 
psia Natural Gasoline between —50°F and +140°F., The intended application of this proce- 
dure is for determining the value of P, (equilibrium vapor pressure at flowing conditions) 
required for determining pressure volume correction factors as specified in API MPMS 
Chapter 11.2.1 and 11.2.2. 

Three primary correlations (commercial HD-5 Propane, Butanes, and Natural Gasoline 
mixtures) were developed from data. To provide continuity, two secondary correlations or 
extrapolations were derived which connect or bridge the three primary correlations. 

Minor variations from the predicted vapor pressure to actual are to be expected in the 
primary equation due to the possible combinations of product compositions which may 
result in the same relative density product. The magnitude of these variations, however, is 
expected to have minimal overall impact on the volume correction factors determined from 
its use. 

Due to the unconstrained range of compositions possible in non-specification mixtures, 
greater variations may occur between an actual sample and the correlation values in the 
extrapolated regions, and accuracy should be checked by the user. 

These bridge correlations may be utilized for settlement purposes if agreed to between 
the contracting parties. However, additional accuracy may be obtained by developing a 
modified correlation as indicated in the procedure, based upon test data. 

API publications may be used by anyone desiring to do so. Every effort has been made 
by the Institute to assure the accuracy and reliability of the data contained in them; how- 
ever, the Institute makes no representation, warranty, or guarantee in connection with this 
publication and hereby expressly disclaims any liability or responsibility for loss or damage 
resulting from its use or for the violation of any federal, state, or municipal regulation with 
which this publication may conflict. 

Suggested revisions are invited and should be submitted to Measurement Coordination, 
Exploration and Production Department, American Petroleum Institute, 1220 L Street, 
N.W., Washington, D.C. 20005. 
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Chapter 11—Physical Properties Data 


ADDENDUM TO SECTION 2, PART 2—COMPRESSIBILITY FACTORS FOR 
HYDROCARBONS, CORRELATION OF VAPOR PRESSURE FOR COMMERCIAL NATURAL 
GAS LIQUIDS 


1 Background 


In custody transfer applications there is a need to convert 
flowing volumes of liquids to equivalent volumes of liquid 
at some set of standard conditions, usually 60°F and the 
greater of one atmosphere pressure or the equilibrium vapor 
pressure of the liquid. The techniques for converting 
densities from flowing conditions to standard conditions 
generally are of the form: 


Poo = XE Xp 
Where: 


P69 = Density at 60°F and the greater of one atmosphere 
pressure or the equilibrium vapor pressure of the 
liquid. 

= Density at flowing conditions. 

Correction factor for temperature effects. 

= Correction factor for pressure effects. 


SILO 


Two methods used for calculation of the F,, term are 
standardized by the American Petroleum Institute (1984 and 
1986). These methods require the equilibrium bubble point 
pressure (P.), the vapor pressure, at the flowing conditions in 
order to calculate the correction factor. Although the vapor 
pressure of the flowing liquid could be measured, it generally 
is not. The vapor pressure can be calculated from composi- 
tional information to the required accuracy using one of sev- 
eral techniques, but the composition is not always measured 
for natural gas liquids, NGLs. Therefore, a correlation for the 
vapor pressure of NGLs based upon normally measured 
properties was desired. 

A correlation for the vapor pressure of NGLs was the cor- 
relation developed by Hillburn (1968) (see Figure 1} which 
requires only the vapor pressure at 100°F. This correlation was 
primarily developed for use with commercial propanes, but it 
has been applied to other NGLs. A major limitation of this cor- 
relation is the fact that it can predict negative absolute vapor 
pressure. Another limitation is that it requires the 100°F vapor 
pressure value. This value is often measured to meet transpor- 
tation regulation requirements, but it is not always measured. 
The correlation seems to work fairly well for propane mix- 
tures, but does not work at all well for natural gasolines. 

The problems and limitations of the Hillburn correlation 
prompted this project to develop an alternative correlation 
based only upon the relative density, a value that is always 
measured. A less empirical correlation form was also desired. 
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Since the vapor pressure curve for a mixture has a shape 
very similar to a pure component vapor pressure curve at low 
to moderate pressures, it was decided to build the correlation 
around one of the simplest pure component vapor pressure 
models that could reproduce experimental vapor pressure 
data reasonably well. Of the several models reviewed, the 
simple Antoine (Reid, et al. 1987) vapor pressure model was 
chosen. It was found that a fairly good correlation of the 
Antoine constants in terms of relative density was possible, 
yielding the final correlation form shown in Figure 2. Table 
] lists the constants required for the correlation in each of its 
applicable ranges. Figure 3 is a plot of the vapor pressures 
that would be predicted by the new correlation at several key 
temperatures in the temperature range of the correlation. 


2 Correlation Development 
2.1 PROPANES 


2.1.1 Product Specifications 


For the purpose of correlation development, product 
specifications from GPA Standard 2140-88, “Liquefied 
Petroleum Product Gas Specifications and Test Methods,” 
given in Table 2, were used. Although the presence of 
olefins is allowed by the specifications, it was assumed that 
no olefins were present since they generally have a much 
higher market value and would be preferentially recovered. 
During the course of the correlation development, it became 


P = Pigg + Q X (Pigg + 20) X Tan(a) 
For t < 100°F 
Q = -1.0 
a = 0.017454 x (36.95865 — 0.1363394 x t 
—0,000927206 x t? - 0.00001409909 x 7) 
For t > 100°F 


Q = 10 

a = 0.0143464 x (t— 100) 
Where: 

P = Vapor pressure at t (psig). 


- 
Il 


Temperature (°F). 
Vapor pressure at 100°F (psig). 


Piao 


Figure 1—Old Correlation for Vapor Pressures 
(Hiliburn, 1968) 
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2 CHAPTER 11—PHYSICAL PROPERTIES DATA 


Table 1—Constants for New Correlation 


Applicable Range Ao Ay Bo B, 
0.490 <G <0510 6.4747 37.083 12038.0 —31296.5 
£0.7023 +1.405 4335.5 +671.0. 
0.510 < G < 0.560 11.5454 1.7490 1378.8 —10396.1 
+1.1657 +1.164 +360.9 +555.9 
0.560 $ G < 0.585 6.4827 10.790 3721.5 -14579.5 
+0.4220 40.735 +201.6 4351.3 
0.585 < G < 0.625 6.5412 10.690 6514.5 -19353.9 
+3.1017 45.201 $518.1 £736.3 
0.625 < G $ 0.676 20.8537 -12.210 6765.6 1894.3 
+3.7280 45,722 +489.0 £750.6 


Note: Constants are in U.S. engineering units for new vapor pressure 
correlation for each relative density range. Constants are given to more 
digits than are significant in order to allow reproduction of results. The 
confidence limits given are 90 percent confidence limits. 


apparent that a propane mix having a relative density of 
0.490 or 0.495 could not meet these specifications. This fact 
was later proven by experimental measurements (Parrish, 
1987). Since several member companies wanted the cortela- 
tion to apply to propane mixes having relative densities of 
0.490 and 0.495, it was agrced that for these cases the L00°F 
vapor pressure specification would be relaxed in order to 
obtain the desired relative density. 


2.1.2 Data Collection 


As part of the correlation project definition, it was agreed 
that the vapor pressure data would be generated from the 
Soave-Redlich-Kwong (SRK) equation of state (Soave, 
1972) and that densities would be generated from the CO- 
STALD correlation (Hankinson, et al., 1979 and Thomson, 
et al., 1982) using parameters from the Phillips Thermo- 
physical Data System (Embry, 1985). It was found that it 
was necessary to modify the values of the acentric factor 
used with the SRK in order to match the 100°F vapor 
pressures to the desired accuracy. The values of the pure 
component parameters used are listed in Table 3. Volume 
percent compositions were converted to molar compositions 
using the Pq values. The binary interaction parameters 
needed for the SRK are shown in Table 4. 

The COSTALD correlation has been shown by Thomson, 
et al. (1982), to give accurate NGL densities, on the order of 
+0.2 percent. A specially modified version of the correlation 
has been developed for LNGs (Hankinson, et al., 1982), but 
was not applicable to this work. 

The compositions of several propane samples which met 
the product specifications at relative densities of 0.490, 
0.495, 0.500, 0.505, and 0.510 were calculated through the 
process of trial and error. It was found that there was enough 
leeway in the specifications that it was possible to generate 
more than one composition that would meet them at each 
relative density. Only the compositions which yielded the 
highest and lowest 100°F vapor pressure, shown in Table 5: 
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were retained for further analysis. The variation between the 
high and low vapor pressure cases for relative densities of 
0.500 and 0.510 was negligible. Therefore, only one case 
was retained for these two relative densities. 

Vapor pressure data was generated for each case over 
the entire temperature range of —50°F to 140°F in 10°F 
increments using the SRK. This provided twenty vapor 
pressure points for each case, for a total of 160 points. 


2.1.3 Data Regression 


For the purpose of data regression it was decided that a 
single set of vapor pressure data would be used at each 
relative density. The “representative” vapor pressure used 
for samples which had both an upper and a lower vapor 
pressure was chosen such that it would yield an equal 
percentage error from each of the true vapor pressures. This 
choice of the “representative” pressure is almost the same as 
would be chosen by the regression. The “representative” 
pressure is calculated as the harmonic average of the high 
and the low vapor pressures. Only the “representative” 
vapor pressure data were used in the data regression. 

Regression of the generated vapor pressure data was 
performed in three steps: 


a. Regression of each sample individually to the equation: 
In(P)=A+B/(t+C). 

b. Regression of A and B as linear functions of relative 

density. 

c. Regression of all samples simultaneously to the correla- 

tion equation. 


By performing the data regression in three steps, a check 
of the appropriateness of the chosen functional forms and a 


[ in (P) 
Where: 


P = Vapor pressure at t (psia). 
t = Temperature (°F). 
G = Relative density (60°F/60°F)*. 
C = Aconstant. For all ranges C = 443.0°F. 
Ao = Values for different correlation ranges are given 
in Table 1. Units are Ln (psia). 
A, = Values for different correlation ranges are given 
in Table 1. Units are Ln (psia). 
Bo = Values for different correlation ranges are given 
in Table 1. Units are Ln (psia)/°F. 
B, = Values for different correlation ranges are given 
in Table 1. Units are Ln (psia)/°F. 


Note: The correlation is limited to 50°F < t < 140°F, and 
0.490 < G < 0.676. 


(Ag + Ay x G) + (By + B, x G)/(+C) 


4 It is assumed that an accurate value of the relative density will be used 
with this correlation. The interested reader should refer to the applicable 
ASTM standards for reliable measurement-methods. 


Figure 2-New Correlation for Vapor Pressures 
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check of the effects of each regression assumption of the 
goodness of fit could be performed. 

Regression of the samples individually yielded fits of the 
data having a maximum error of 1.0 percent and an average 
absolute error of 0.3 percent over the entire temperature 
range. A residual plot did not show any significant system- 
atic variations, so the correlation form was accepted without 
any further investigation. 

At this point, several different values of C were tried to see 
if the maximum and average errors could be reduced. There 
was no significant reduction of error possible, so the value of 
C was “frozen” at the reported value of 443°F for all further 
regression analyses. 

Plots of the A and B values determined by the individual 
regressions as functions of relative density showed an almost 
linear dependence. Residual plots from the linear regressions 
of the A and B values showed some systematic variation, but 
addition of higher order relative density terms neither 
improved the fit nor removed the systematic deviations. 
Therefore, the linear dependence was chosen. 


The previous sets of analyses yielded the correlation form 
given in Figure 2 which was then used for regression of all 
the propane data. For the data included in the analysis, the 
maximum error was 3.9 percent and the average absolute 
error was 0.8 percent. When checked against all of the data 
rather than just the “representative” data, the correlation had 
the maximum error of 7.8 percent, an average absolute error 
of 1.8 percent and a bias of 0.2 percent. The ability of the 
correlation to reproduce the raw propanes data at —50°F, 0°F, 
60°F, 00°F, and 140°F is shown in Figure 4, Figure 5, 
Figure 6, Figure 7, and Figure 8 respectively. 


2.1.4 Correlation Comparison 


Comparison of the new correlation to the Hillburn corre- 
lation for the propanes shows that the new correlation is 
superior. Table 6 gives a summary comparison of the new 
correlation and the old Hillburn correlation at selected points 
for each of the samples. At most temperatures the new corre- 
lation is superior. 


Table 2—Liquefied Petroleum Gas Specifications (GPA Standard 2140-88) 


Product Commercial Commercial Commercial B-P 
Characteristics Propane Butane Mixtures Propane HD-5 Test Method 
Predominantly Not less than 90 
mixtures of liquid volume percent 
Predominantly Predominantly butanes and/or propane; not more than 
propane and/ butanes and/ butytenes with propane five liquid volume ASTM 
Composition or propylene or butylenes and/or propytene percent propylene D-2163-82 
Vapor pressure 208 70 208 208 ASTM 
at 100°F, psig, D-1267-84 
maximum 
Temperature at 95 -37 36 36 -37 ASTM 
percent evaporation, D-1837-86 
°F maximum 
Butane or heavier, 2.5 a - 2.5 ASTM 
liquid volume D-2163-82 
percent, maximum 
Pentane or heavier, - 2.0 2.0 - ASTM 
liquid volume D-2163-82 
percent, maximum 
Residue on 0.05 ml - oa 0.05 ml ASTM 
evaporation of 100 D-2158-85 
milliliters, maximum 
Oil stain observation Pass* - Pass* ASTM 
D-2158-85 
Corrosion, copper No. | No. | No. | No. | ASTM 
strip, maximum D-1838-84 
Total sulfur, ppmw 185 140 140 123 ASTM 
D-2784-80 
Moisture content Pass - - Pass GPA Propane 
Dryness Test 
(Cobalt Bromide) 
or D-2713-86 
Free water content - None - None - 


* An acceptable product shall not yield a persistent oil ring when 0.3 milliliters of solvent residue mixture is added to a filter paper in 0.1 increments and 


examined in daylight after two minutes as described in ASTM D-2158. 
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Table 3—Correlation Constants for Propanes and Butanes 


Component T. (PF) P. (psia) ® used with SRK SRK V#*(cm3/mol} Péo (b/ft?) 
Ethane 90.09 707.80 0.0983 0.0983 145.80 22.23 
Propane 206.01 616.30 0.1667 0.1532 200.10 31.66 
n-Butane 305.65 550.70 0.2008 0.2008 254.40 36.45 
i-Butane 274.98 529.10 0.1867 0.1825 256.80 35.12 
n-Pentane 385.70 488.60 0.2522 0.2522 311.30 39.35 


Note: These constants were used when developing the vapor pressure correlations for commercial propanes and commercial butanes. The constants were taken 
from the Phillips Thermophysical Data System. The acentric factors (a) for vapor pressure were adjusted to give a better fit of the 100°F vapor pressures. 


Table 4—SRK Interaction Parameters for Propanes and Butanes 


Component j 
Component i Ethane Propane n-Butane ji-Butane n-Pentane 
Ethane 0.0 0.002 0.005 0.0 0.007 
Propane 0.002 0.0 0.0 0.0 0.003 
n-Butane 0.005 0.0 0.0 0.0 0.0 
i-Butane 0.0 0.0 0.0 0.0 0.0 
n-Pentane 0.007 0.003 0.0 0.0 0.0 


Note: These parameters were used in the SRK correlation for calculation of the propanes and butanes vapor pressures. 


Table 5—Compositions and Relative Densities of Propane Samples 


Volume Percent of Component 


Case Name Relative Density (60°F/60°F) Ethane Propane n-Butane i-Butane n-Pentane 
0.490 G Hi V.P. 0.49000 20.00 77.50 2.45 0.00 0.05 
0.490 G Lo VP. 0.49000 17.70 82.30 0.00 0.00 0.00 
0.495 G Hi V.P. 0.49500 15.20 82.30 2.45 0.00 0.05 
0.495 G Lo V.P. 0.49500 12.73 Bh27 0.00 0.00 0.00 
0.500 G HD-5 0.50036 7.35 92.65 0.00 0.00 0.00 
0.505 G Hi V.P. 0.50500 4.81 92.69 0.45 0.00 0.05 
0.505 G Lo VP. 0.50500 2.36 97.64 0.00 0.00 0.00 
0.510 G HD-5 0.50950 0.00 97.50 2.45 0.00 0.05 


Note: These are actual compositions and calculated relative densities used for development of the correlation in the commercial propanes relative density range. 
Only the cases having a relative density of 0.500 or greater met the 100°F vapor pressure requirement for commercial propanes. 
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Table 6é—Comparison of Vapor Pressure Correlations for Commercial Propanes 


New Correlation 


Case Name 


Relative 
Density 
(60°F/60°F) 


“Representative” 


0.490 G 


0.490 G Hi V.P. 


0.490 Lo V.P. 


“Representative” 


0.495 


0.495 Hi VP. 


0.495 Lo V.P. 


0.500 G HD-S 


“Representative” 


0.505 


0.505 Hi V.P. 


0.505 Lo V.P. 


0.510G HD-5 


Average Percent Error 


0.4900 


0.4900 


0.4900 


0.4950 


0.4950 


0.4950 


0.5003 


0.5050 


0.5050 


0.5095 


Average Absolute Percent Error 


Maximum Percent Error 
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Temperature 
(°F) 


SRK Vapor 
Pressure 
(psia) 


26.87 
68.92 
168.58 
275.84 
423.16 


26.05 
67.27 
165.63 
272.08 
418.97 


27.74 
70.66 
171.63 
279,71 
427.43 


22.92 
60.49 
151.94 
252.29 
392.21 


22.05 
58.72 
148.78 
248.21 
387.57 


23.86 
62.36 
155.24 
256.61 
396.97 


17.73 
49,52 
130.76 
222.75 
354.10 


14.54 
42.72 
117.25 
203.38 
327.99 


15.47 

44.54 
120.41 
207.40 
332.54 


13.72 
41.04 
114.25 
199.51 
323.56 


11.59 
36.36 
104.54 
185.14 
303.35 


Vapor 


Pressure 


(psia) 


27.27 
70.30 
170.83 
276.87 
419.97 


27.27 
70.30 
170.83 
276.87 
419.97 


27.27 
70.30 
170.83 
276.87 
419.97 


22.05 
59.44 
150.65 
249,83 
386.52 


22.05 
59.44 
150.65 
249.83 
386.52 


22.05 
59.44 
150.65 
249.83 
386.52 


17.58 
49.72 
131.80 
223.96 
353.86 


14.41 
42.49 
117.17 
203.42 
327.41 


14.41 
42.49 
117,17 
203.42 
327.41 


14.41 
42.49 
117.17 
203.42 
327.41 


11.89 
36.54 
0.09 
185.45 
303.84 


Percent 


Error 


Old Correlation 


Vapor 
Pressure 
(psia) 


11.64 
64.30 
160.00 
275.84 
437.18 


11.42 
63.37 
157.78 
272.08 
431.25 


11.88 
65.26 
162.27 
279.71 
443.27 


10,22 
58.47 
146.15 
252.29 
400.11 


9.98 
57.46 
143.75 
248.21 
393.69 


10.48 
59.52 
148.63 
256.51 
406,75 


8.44 
51.16 
128.78 
222.75 
353.61 


7.28 
46.36 
117.39 
203.38 
323.13 


7.52 
47.36 
119.75 
207.40 
329.46 


7.04 
45.41 
115.12 
199.51 
317.05 


6.18 
41.85 
-2.04 

185.14 
294.43 


Percent 
Error 
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2.2 BUTANES 


2.2.1 Product Specifications 


Specifications given by the Gas Processors Association 
for butanes are shown in Table 2. For the purposes of corre- 
lation development, it was decided to specify two classes of 
commercial butane: 


a. A normal butane product having at least 96.5 volume 
percent n-butane and containing no more than 3.5 volume 
percent i-butane. 

b. An i-butanc product having at least 95.5 volume percent 
i-butane and containing no more than 4.5 volume percent 
normal butane or 2.5 volume percent propane. 


All correlation development was based around these two 
butane product specifications. 


2.2.2 Data Collection 


The same general data collection techniques as were used 
for the propanes were used for the butanes. The parameters 
used in the correlation development are shown in Table 3 and 
Table 4. 

A different approach for generating the compositions was 
used for the butane products than for the propane products. 
Five different compositions were used for each product: 


a. A low vapor pressure case containing the maximum 
allowable pentanes. 

b. A low vapor pressure case containing the minimum 
amount of the primary component and no lights. 

c. A pure component case. 


d. A high vapor pressure case containing the minimum 
amount of the primary component and no heavies. 

e. A typical commercial product composition supplied by 
one of the member companies. 


The compositions of generated cases are given in Table 7. 
Vapor pressure data for each case was generated at 10° incre- 
ments over the temperature range of ~5O°F to 140°F, This 
provided 200 points to be included in the regression. 


2.2.3. Data Regression 


The same data regression procedure was used for fitting 
the butanes data, except no “representative” data was gener- 
ated. All of the data generated was used when fitting the 
correlation. When each case was fitted individually, the av- 
erage absolute error was 0.9 percent, and the maximum error 
was 2.7 percent. The fit of all data to the butanes correlation 
gave an average absolute error of |.1 percent and a maxi- 
mum error of 4.5 percent. The ability of the correlation to re- 
produce the butanes data at -50°F, O°F, 60°F, 100°F, and 
140°F is shown in Figure 9, Figure 10, Figure 11, Figure 12, 
and Figure 13 respectively. 


2.2.4 Correlation Comparison 


Comparison of the new correlation to the Hillburn 
correlation for the butanes shown that the new correlation is 
superior. Table 8 gives a summary comparison of the new, 
proposed correlation and the old, Hillburn correlation at 
selected points for each of the samples. At most temperatures 
the new correlation is superior. 


Table 7—Compositions and Relative Densities of Butane Samples 


Volume Percent of Component 


Case Name Relative Density (60°F/60°F) Ethane 
i-C4 Lo V.P-1 0.5643 0.00 
i-C4 Lo V.P-2 0.5639 0.00 
i-C4 Pure 0.5628 0.00 
i-C4 Hi VP. 0.5616 0.00 
i-C4 Koch 0.5624 0.00 
n-C, Lo V.P-1 0.5854 0.00 
n-C, Lo V.P-2 0.5852 0.00 
n-Cy Pure 0.5844 0.00 
n-C, Hi V.P. 0.5837 0.00 
n-Cy Koch 0.5843 0.00 


Propane n-Butane i-Butane n-Pentane 
0.00 0.00 98.00 2.00 
0.00 4.50 95.50 9.00 
0.00 0.00 100.00 0.00 
2.50 0.00 97.50 0.00 
1.70 1.50 96.80 0.00 
0.00 98.00 0.00 2.00 
0.00 98.50 0.00 1.50 
0.00 100.00 0.00 0.00 
0.00 96.50 3.50 0.00 
0.00 97.50 2.00 0.50 


Note: These are actual compositions and calculated relative densities used for development of the correlation in the commercial butanes’ relative density range. 
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Table 8—Comparison of Vapor Pressure Correlations for Commercial Butanes 


Relative 
Density 
Case Name (60°F/60°F) 


i-C, Hi V.P. 0.5616 


i-C4 Koch 0.5624 


i-C4 Pure 0.5628 


i-C4 0.5639 


iH-Cy 0.5643 


n-Cy Hi VP. 0.5837 


n-Cy Koch 0.5843 


n-Cy Pure 0.5844 


n-C4 Lo V.P-2 0.5852 


n-Cy Lo V.P-1 0.5854 


Average Percent Error 
Average Absolute Percent Error 


Maximum Percent Error 
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Temperature 
(CF) 


-50 


SRK Vapor 
Pressure 
(psia) 


3.26 
12.01 
39.77 
75.66 

131.43 


3.16 
11.71 
38.99 
74.35 

129.41 


3.01 
11.28 
37.90 
75.28 

126.78 


2.95 
11.09 
37.34 
71.61 

125.19 


2.96 
11.10 
37.30 
71.46 

124.81 


1.77 
7.25 
26.32 
52.44 
94.58 


1.74 
7.16 
26.06 
51.97 
93.81 


1.73 
7.10 
25,91 
51.74 
93.46 


1.7] 
7.03 
25.65 
51.23 
92.57 


1.70 
7.00 
25.56 
51.07 
92.28 


New Correlation 


Vapor 
Pressure Percent 
(psia) Error 
3.25 0.40 
11.71 2.50 
38,98 1.99 
74.97 0.91 
131.82 0.30 
3.18 -0.61 
11.50 1.76 
38.41 1.48 
74.01 0.45 
130.33 0.71 
3.15 4.52 
11.40 -1,09 
38.14 -0.62 
73.54 ~1,32 
129.59 -2.22 
3.06 -3.60 
11.13 0.35 
37.38 -0.11 
72,25 -0.90 
127.58 -1.9) 
3.02 2.17 
11.03 0.63 
37.11 0.5] 
71.79 0.46 
126.85 -1.64 
1.82 -2.56 
7.18 0.94 
26.07 0.94 
52.57 0.25 
96.28 -1.79 
1.79 -2.70 
7.09 1.02 
25.79 1.04 
$2.07 0.19 
95.46 -1.76 
1.78 3.02 
707 0.41 
25.74 0.65 
51.98 0.47 
95.32 —1.99 
1.75 —2.05 
6.95 1.18 
25.37 1.08 
51.32 0.18 
94,25 ~1.81 
1.74 —2.12 
6.92 1.19 
25,28 1.10 
31.16 0.17 
93.97 -1.84 
0.01 
1.10 
4.52 


Old Correlation 


Vapor 
Pressure 
(psia) 


0.42 
14.74 
42,30 
75,66 
122.12 


-0.50 
14.42 
41.53 
74.35 
120,06 


-0.61 
13.98 
40.49 
72.58 
117.27 


0.67 
13.74 
39.92 
71,61 
115.75 


—0.68 
13.70 
39.83 
71.46 
115.51 


-1.82 

8.99 
28.65 
52,44 
85.58 


-1.85 

8.88 
28.37 
51.97 
84.84 


-1.87 

8.82 
28.24 
31.74 
84.47 


-1.90 

8.69 
27.93 
51.23 
83.67 


-1.91 

8.65 
27.84 
$1.07 
83.42 


Percent 
Error 


112.98 
22.74 
6.36 
0.00 
7,08 


115.89 
—23.12 
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2.3. NATURAL GASOLINES 


2.3.1 Product Specifications 


According to the GPA Standard 3132-84, “Natural Gas- 
oline Specifications and Test Methods,” natural gasoline is 
defined for commercial purposes by the specifications in 
Table 9 as follows: 


Table 9—Natural Gasoline Specifications 
(GPA Standard 3132-84) 


Product Characteristic Specification Test Method 

Reid Vapor Pressure 10-34 pounds ASTM D-323 

Percentage Evaporated 25-85 ASTM D-216 

at 140°F 

Percentage Evaporated not less than 90 ASTM D-216 

at 275°F 

End Point not more than ASTM D-216 

375°F 

Corrosion not more than ASTM D-130 
classification 1 (modified) 

Color not less than plus ASTM D-156 

25 (Saybolt) 
Reactive Sulfur Negative, “sweet” GPA 1138 


There are also 24 different grades of natural gasoline defined 
by six ranges of Reid Vapor Pressure and four ranges of per- 
centage evaporated at 140°F as shown in Table 10: 


Table 10—Grades of Natural Gasoline 


Liquid Volume Percentage Evaporated at 140°F 


Reid Vapor 225&<40 240&<55 255&<70 270k <85 
Pressure (psia) Grade Grade Grade Grade 
230 & <34 34-25 34-40 34-55 34-70 
>26 & = 30 30-25 30-40 30-55 30-70 
>22 & $26 26-25 26-40 26-55 26-70 
>18 & <22 22-25 22-40 22-55 22-70 
>14&<18 18-25 18-40 18-55 18-70 
Z10&< 14 14-25 14-40 14-55 14-70 


Note: The various grades of natural gasoline are as specified by the Gas 
Processors Association. 
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The large variance in composition, vapor pressure, and 
relative density allowed by the many grades of natural gaso- 
line makes development of a single vapor pressure correla- 
tion for natural gasoline very difficult. At a constant relative 
density, it is possible to adjust the composition to yield as 
much as a 40 percent variation in the 100°F vapor pressure. 
This means that a vapor pressure correlation based upon 
relative density alone cannot be expected to predict the 
vapor pressure of a single sample better than + 20 percent. 


2.3.2 Data Collection 


Based upon the specifications for natural gasolines, it was 
determined that many more components were necessary for 
characterizing natural gasoline than were used for charac- 
terizing the propanes and butanes. The pure component 
parameters used for generating the necessary data are listed 
in Table 11. The corresponding SRK binary interaction 
parameters are listed in Table 12. Again, all of this data was 
taken from the Phillips Thermophysical Data System. 


The large variation of properties allowable for natural gas- 
olines made generation of compositions that meet the speci- 
fications an intractable task. Therefore, representative 
natural gasoline compositions were supplied by member 
companies (see Table 13) and were used for data generation. 
Vapor pressure data was then generated at 10°F increments 
over the temperature range —50°F to 140°F using the SRK 
equation of state. 


2.3.3 Data Regression 


The stepwise data regression procedure outlined for the 
propanes was also used for the natural gasolines. When each 
case was regressed individually, the average absolute error 
was 1.2 percent, and the maximum error was 3.9 percent. It 
might have been possible to reduce these errors by adjust- 
ment of the C constant, but that was not attempted. 


A plot of the regressed A and B constants showed a slight 
linear dependence upon specific gravity, but there was 
considerable scatter. Simultaneous regression of all of the 
data gave very large statistical confidence limits for the 
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regression parameters. The average absolute error was 14.4 
percent, and the maximum error was over 54 percent. These 
kind of errors were expected since that kind of variation was 
present in the data. The ability of the correlation to reproduce 
the natural gasolines data at —50°F, O°F, 60°F, 100°F, and 
140°F is shown in Figure 14, Figure 15, Figure 16, Figure 
17, and Figure 18 respectively. 

When the data were regressed simultaneously, but with a 
different A constant for each case, the By and B, were 
almost identical to those obtained from the simultaneous 
regression. However, the average absolute error was 
reduced to 3.6 percent and the maximum error was reduced 
to 16.8 percent. This dramatic reduction in error gives 
credence to the ideas outlined in “Improved Accuracy.” 


2.3.4 Correlation Comparison 


Comparison of the new correlation to the Hillburn corre- 
lation for the natural gasolines shown that the new correla- 
tion is superior. Table 14 gives a summary comparison of 
the new, proposed correlation and old, Hillburn correlations 
at selected points for each of the samples. At most tempera- 
tures the new correlation is superior. The Hillburn correla- 
tion is especially poor for these samples in that it predicts 
negative absolute pressures for many samples at the lower 
temperatures. Although the new correlation is not nearly as 
good for these samples as it was for the propanes and 
butanes, it is still considerably better than the Hillburn 
correlation. 


Table 11—Correlation Constants for Natural Gasolines 


Component T, OF) P. (psia) «w used with SRK SRK V#(cm*/mol) Dgo (lb/ft?) 
n-Butane 305.65 550.70 0.2008 0.2008 254.40 36.45 
i-Pentane 369.10 490.40 0.2400 0.2400 309.60 38.96 
n-Pentane 385.70 488.60 0.2522 0.2522 311.30 39.35 
n-Hexane 453.70 436.90 0.3007 0.3007 368.20 4141 
n-Heptane 512.80 396.80 0.3507 0.3507 430.40 42.92 
n-Octane 564,22 360.60 0.3998 0.3998 490.40 44.08 
n-Nonane 610.68 332.00 0.4478 0.4478 552.90 45.01 
n-Decane 652.10 304.00 0.4916 0.4916 619.20 45.79 
Methylcyclohexane 570.27 503.50 0.237] 0.2371 370.90 48.27 


Note: These constants were used when developing the vapor pressure correlations for natural gasolines. The constants were taken from the Phillips Thermo- 
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physical Data System. 


Table 12—SRK Interaction Parameters for Natural Gasolines 


Component j 
Methylcyclo- 

Component i n-Butane i-Pentane n-Pentane n-Hexane n-Heptane n-Octane n-Nonane n-Decane hexane 
n-Butane 0.0 0.0 0.0 0.003 0.005 0.007 0.010 0.010 0.0 
i-Pentane 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 
n-Pentane 0.0 0.0 0.0 0.0 0.003 0.005 0.009 0.009 0.0 
n-Hexane 0,003 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 
n-Heptane 0.005 0.0 0.003 0.0 0.0 0.0 0.0 0.0 0.0 
n-Octane 0.007 0.0 0.005 0.0 0.0 0.0 0.0 0.0 00 
n-Nonane 0.010 0.0 0.009 0.0 0.0 0.0 0.0 0.0 0.0 
n-Decane 0.010 0.0 0.009 0.0 0.0 0.0 0.0 0.0 00 
Methylcyclo- 

hexane 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 


Note: These parameters were used in the SRK correlation for calculation of the vapor pressures of natural gasolines. 
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3 Extension and Improvement of 
Correlation 


3.1. INTERMEDIATE COMPOSITIONS 


Extensions of the correlations developed for propanes, 
butanes, and natural gasolines to intermediate compositions 
or blends of these products were developed. It was assumed 
that the same correlation form could be used for the blends 


3.2 IMPROVED ACCURACY 


Higher accuracy estimates of vapor pressures at tempera- 
tures other than 100°F can be obtained from the correlation 
by use of an experimental 100°F vapor pressure point. The 
improvement is most significant for the natural gasolines 
where the variation of 100°F vapor pressure of products 
having the same relative density was most significant. The 
modified form of the equation is shown in Figure 19. 


as for the original products. This is a relatively safe assump- 
tion, since the same correlation form worked adequately for 
all products considered. It was further assumed that the 
blends would have a relative density dependence that was 
intermediate to the bracketing products. This assumption 
would only be valid when the blend is composed solely of 
the bracketing products. The problems encountered with 
correlating natural gasolines showed that strict composition 
restrictions must be placed on the products in order for a 
satisfactory correlation to be developed. 

Based upon the preceding assumptions, the intermediate 
composition correlations were developed solely from the 
correlations of the bracketing products, and no new infor- 
mation was used to check the correlations at the intermedi- 
ate compositions. The values of the correlation constants 
were determined such that the pressures predicted by the 
interpolating correlation at its limits would exactly match 
the predicted pressures of the bracketing correlations at the 
abutting limits. 


In(P) = In(P 199) + (Bo + By x G) x (100.0 — t)(D x (t+ C)) 
Where: 


P = Vapor pressure at t (psia). 
Pigg = Vapor pressure at 100°F. 
t = Temperature (°F). 
G = Relative density (60°F/60°F). 
C = Aconstant. For all ranges C = 443.0°F. 
D = Aconstant. For all ranges D = 543.0°F. 
Bo = Values for different correlation ranges are given 
in Table 1. Units are Ln (psia)/°F. 
B, = Values for different correlation ranges are given 
in Table 1. Units are Ln (psia)/°F. 


Note: The correlation is limited to -SO°F < ts 140°F, and 
0.490 < G < 0.676. 


Figure 19—Modified Correlation for Vapor Pressures 


Table 13—Compositions and Relative Densities of Natural Gasolines 


Valume Percent of Component 


Relative Methyl- 

Case Density n- i- n- n- n- n- n- n- cyclo- 

Name (60°F/60°F) = Butane Pentane Pentane Hexane Heptane Octane Nonane Decane hexane 
MAPCO 0.6672 6.00 44.00 0.00 25.00 6,00 1.50 0.00 0.00 17.50 
Max V.P. 
MAPCO 0.6693 0.00 50.00 0.00 25.00 6.00 1.50 0.00 0.00 17.50 
Min V.P. 
MAPCO 0.6601 3.53 25.43 30.79 21.20 4.58 1.14 0.00 0.00 13.30 
#12 
Cities 0.6760 0.83 24.29 23.42 19.72 7,57 1.89 0.00 0.00 22.07 
Service 
Warren 1-Cs 0.6246 5.00 95.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
Max V.P. 
Warren i-Cs 0.6269 0.00 97.50 2.50 0.00 0.00 0.00 0.00 0.00 0.00 
Min V.P. 
Warren n-Cs 0.6306 0.00 5.00 95.00 0.00 0.00 0.00 0.00 0.00 0.00 
Max V.P. 
Warren n-Cs 0.6327 0.00 0.00 95.00 5.00 0.00 0.00 0.00 0.00 0.00 
Min VP. 
Koch M 
Grade 0.6656 4.63 23.09 25.89 23.20 5.57 1,39 0.00 0,00 16.24 
Koch 
Low V.P. 0.6588 17.88 18.63 19.41 22.10 5.30 1.33 0.00 0.00 15.46 


Note: These are the actual compositions and calculated relative densities for use for development of the correlation in the natural gasolines’ relative density range. 
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Table 14—-Comparison of Vapor Pressure Correlations for Natural Gasolines 


Relative 
Density 
Case Name (60°F/60°F) 
Warren 1-Cs 0.6246 
Max V.P. 
Warren i-Cs 0.6269 
Min V.P. 
Warren n-Cs5 0.6306 
Max V.P. 
Warren n-Cs 0.6327 
Min V.P. 
Koch Low V.P. 0.6588 
MAPCO #12 0.6601 
Koch M Grade 0.665672 
MAPCO Max V.P. 0.6672 
MAPCO Min V.P. 0.6693 
Cities Service 0.6760 


Average Percent Error 
Average Absolute Percent Error 


Maximum Percent Error 
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Temperature 
(°F) 


SRK Vapor 
Pressure 


(psia) 


0.47 
2.13 
9.94 
21.62 
41.90 


0.38 
1.98 
8.86 
19.60 
38.51 


0.27 
1.46 
6.88 
15.66 
31.48 


0.25 


New Correlation 


Vapor 


Pressure 


(psia) 


0.38 
[.87 
8.4] 
19.05 
38.57 


0.37 
1.84 
8.25 
18.67 
37.78 


0.36 
1.78 
7.99 
18.08 
36.55 


0.35 
1.75 
7.85 
17.75 
35.87 


0.29 
1.43 
6.30 
14.14 
28.39 


0.29 
1.41 
6.23 
13.98 
28.06 


0.28 
1.35 
5.95 
13.32 
26.71 


0.27 
1.33 
5.87 
13.14 
26.33 


0.27 
1.3 
5.77 
12.90 
25.84 


0.26 
1.24 
5.45 
12.17 
24.33 


Percent 
Error 


19.84 
19.12 
15.41 
11.90 

7.96 


3.41 
7.36 
6.86 
4.75 
1.89 


35.84 
—22.25 
-16.21 
-15.42 
—16.08 


—42.35 
~27.11 
19.87 
-18.57 
18.85 


41.93 
36.50 
29.3) 
24,31 
19.21 


-7.14 
-3.19 
-2.43 
-3.47 
-5.25 


—1.56 
0.50 
0.04 

-1.44 

-3.$3 


13.94 
13.19 
10.34 
7.61 
4.52 


24,23 
~13.80 
8.94 
-8.30 
~8.83 


-32.46 
—20.22 
—13.88 
-12.3] 
—12.05 


-1.59 
14.16 
42.35 
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Old Correlation 


Vapor 
Pressure 
(psia) 
3.68 

1.36 


Percent 
Error 


883.69 
41.09 
-5.89 

0.00 
11.53 


1092.34 
56.52 
5.46 

0.00 
11.99 


1623.90 
107.74 
—2.09 
0.00 
12.04 


1738.63 
120.61 
0.97 


147.25 


1631.33 
154.45 
6.99 
0.00 
8.08 


1392.48 
130.58 
5.66 
0.00 
7.80 


2058.82 
190.34 
8.98 
0.00 
8.21 


2327.91 
226.57 


188.71 
2327.91 
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APPENDIX A—SAMPLE IMPLEMENTATION 


The following FORTRAN subroutine is a sample implementation of the method outlined in this report. This report imple- 
mentation incorporates several proposed rounding rules that should allow four decimal places of accuracy for any metering 
correction to be obtained. 


SUBROUTINE VAPRS (T, RDEN, PE, VP) 


KKEKKKKKKKKKEKRKKRK KEK KEKE KEK KEK KEKE KK KKEKKEKEKKEKKKKEKKEKKEKKEKEKKEKKKKKKKK KK KK 


THIS SUBROUTINE ESTIMATES THE VAPOR PRESSURE AT METER 
TEMPERATURE USING RELATIVE DENSITY INPUT AND GPA TP-15 


THIS CORRELATION IS VALID FOR 

A RELATIVE DENSITY RANGE OF 0.490 TO 0.676 INCLUSIVE AND 
A TEMPERATURE RANGE OF -50°F TO 140°F INCLUSIVE. 

INPUT DATA BEYOND THESE SPECIFIED RANGES WILL RESULT 

IN PE1 AND VPA BEING SET TO 99999. 


T ---- OBSERVED TEMPERATURE, DEGREES °F 

RDEN - RELATIVE DENSITY 

PE1 -~ CALCULATED VAPOR PRESSURE AT TEMPERATURE T (PSIA) 
VPA -- CALCULATED VAPOR PRESSURE AT 100°F (PSTIA) 


qAeaqqnqnnnnannnnnnnn 


DOUBLE PRECISION T, RDEN, PE, VP 
DOUBLE PRECISION A, B, C, D, NVP, NPE 
INTEGER IPE, IVP, K, L, M 


CHECK THAT INPUT DATA IS WITHIN CORRELATION LIMITS 
IF NOT KICK OUT OF CALCULATIONS 


AANA 


IF (RDEN .LT. 0.490D0} GO TO 1111 
IF (RDEN .GT. 0.676D0} GO TO 1111 
IF (T .LT. —50.0D0) GO TO 1111 
IF (T .GT. 140.0D0) GO TO 1111 


VALID DATA CASE, SO DETERMINE WHICH SET OF CONSTANTS TO USE 


QA 


IF (RDEN .GT. 0.510D0) GO TO 10 
Cc CONSTANTS FOR 0.490 .LE. RDEN .LE. 0.510 


(CONTINUED ON NEXT PAGE) 


21 


Copyright por American Petroleum Institute 
Fri Jan 11 21:48:56 2002 


API MPMS*L1.2.2 ADDENDUM 94 MM 0732290 O542173 T54 mes 


22 CHAPTER 11—PHYSICAL PROPERTIES DATA 
Cc 
A = ~-6.4747D0 
B = 3.7083D1 
C = 1.20380D4 
D = -3.12965D4 
GO TO 50 
Cc 
10 IF (RDEN .GE. 0.560D0) GO TO 20 
Cc 
Cc CONSTANTS FOR 0.510 .LT. RDEN .LT. 0.560 
Cc 
A = 1.15454D1 
B = 1.7490D0 
C.= 1.3788D3 
D -1.03961D4 
GO TO 50 
Cc 
20 IF (RDEN .GT. 0.585D0) GO TO 30 
Cc 
C CONSTANTS FOR 0.560 .LE. RDEN .LE. 0.585 
Cc 
A = 6.4837D0 
B = 1.0790D1 
Co t= 3.7215D3 
D = -1.45795D4 
GO TO 50 
Cc 
30 IF (RDEN .GE. 0.625D0) GO TO 40 
1S 
C CONSTANTS FOR 0.585 .LT. RDEN .LT. 0.625 
Cc 
A = 6.5412D0 
B= 1.0690D1 
Cc = 6.5145D3 
D -1.93539D4 
GO TO 50 
Cc 
40 CONTINUE 
Cc 
Cc CONSTANTS FOR 0.625 .LE. RDEN .LE. 0.676 
Cc 
A = 2.085371D1 
B = -1.2210D1 
C = -6.7656D3 
D = 1.8943D3 
GO TO 50 


(CONTINUED ON NEXT PAGE) 
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Cc 
50 CONTINUE 
Cc 
G CALCULATE THE PROPERLY ROUNDED EXPONENT VALUES FOR PE AND VP 
C 
K = (C + D*RDEN) / (443.0D0 + T) * 1.0D5 
M = B*RDEN * 1.0D5 
L = (K+M+5)/10 
NPE = L 
NPE = A + NPE/1.0D4 
Cc 
K = (C + D*RDEN) / 553.0D0 * 1.0D5 
M = B*RDEN * 1.0D5 
G = (K+M+5) / 10 
NVP = L 
NVP = A + NVP / 1.0D4 
Cc 
Cc CALCULATE VAPOR PRESSURE ROUNDED TO 1 DECIMAL 
Cc 
PE = DEXP (NPE) 
IPE = PE*1.0D1 + 0.5D0 
PE = IPE/1.0D1 
Cc 
VP = DEXP(NVP) 
IVP = VP*1.0Di + 0.5D0 
VP = IVP / 1.0D1 
RETURN 
Cc 
Cc ****X*XTNPUT RELATIVE DENSITY OR TEMPERATURE IS OUT OF RANGE 
Cc **k*k*x*PR AND VP IS SET TO 999999 TO FLAG ANY ERRONEOUS ANSWER 
Cc **k*k*XEROM POSSIBLY BEING USED IN FURTHER CALCULATIONS 
Cc 


1111 CONTINUE 
PE = 999999.0D0 
VP = 999999.0D0 
RETURN 
END 
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